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Abstract 

The spectra of the QSOs PKS2126-158 {z em = 3.27) and Q0055-269 (z em = 3.66) have 
been observed at high resolution (2.2 — 2.7 X 10 4 ). Significant clustering, with £ ~ 1 at 
Av = 100 km s _1 , is detected in both cases for lines with log Njji > 13.8. A similar result 
is obtained analysing the data available in the literature for Q0014+81. An anticlustering 
signal is also detected on scales Av ~ 600 km s _1 . Two voids of size ~ 2000 km s _1 are 
found in the spectrum of Q0055-269 with a probability of 2 X 10~ 4 . 



1 Clustering properties of the Ly a lines 

No clustering in the velocity space has been detected so far for Ly a lines on scales 300 < Av < 
30000 km s -1 ([11], [12], [15]). Their spatial distribution is different from that of the metal-line 
systems selected by means of the CIV doublet, that are known to cluster on scales at least up 
to 600 km s -1 ([13]). Preliminary results at high resolution seem to indicate weak clustering 
on smaller scales (Av = 50 - 300 km s _1 , [14], [1]). 

Although clustering is expected to decrease with redshift, from an empirical point of view it 
is attractive to push the observations towards high redshifts: the larger the density of lines, the 
better the chances of detecting a clustering signal. For this reason we have chosen two relatively 
high-z QSOs, PKS216-158 (z em = 3.27) and Q0055-269 (z em = 3.66), to be observed at ESO 
La Silla with the NTT telescope and the EMMI instrument in the echelle mode ([2], [5] and 
unpublished observations). The weighted mean of the flux-calibrated spectra have a resolution 
R ~ 27000 and R ~ 22000 for PKS2126 and Q0055, respectively. The signal-to-noise ratio per 
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Figure 1: Two-point correlation function of the Lyman— a lines: PKS2126-158 (dot-dashed 
line), Q0055-269(dotted line), Q0014+81 (dashed line) and weighted mean of the three (con- 
tinuous line). For each histogram the horizontal continuous line defines the zero level. The 
distance between two tick marks on the y-axis corresponds to 0.5 units. The continuous curve 
superimposed to the lowest histogram is the model described in Sec. 2 computed from eq. (2) 
with 7 = 1.77, a = 150 km s -1 , r c i = 110 kpc and r Q = 280 kpc at z = 3. 



resolution element at the continuum level ranges from S/N ~ 10 to 40 in the region between Ly a 
and Ly^. Column densities and Doppler widths of the Ly a lines and metal- line systems have 
been derived by a fitting procedure ([6]). The statistical distribution of the Doppler parameter 
for the Ly a lines is peaked at b ~ 23 km s -1 . The column density distribution is described by a 
power-law with a break or cutoff at log Nm ~ 14.5. A featureless distribution is rejected with 
a probability of > 99.9%. 

We have computed the two-point correlation function in the velocity space for PKS2126-158, 
Q0055-269 and Q0014+81 ([10]) as: 



£(v,Av) 



N obs (v,Av) 
N exp (v,Av) 



(1) 



where N {, s is the observed number of line pairs with velocity separations between v and v + Av 
and N exp is the number of pairs expected from a random distribution in redshift. In our line 
sample N exp has been obtained averaging 1000 numerical simulations generated according to 



the cosmological distribution oc (1 + z) 7 , derived from the maximum likelihood analysis of the 
real line sample in the appropriate interval of column density and redshift. In this way it is 
possible to correct for incomplete wavelength coverage due to gaps in the spectra or occultation 
of weak lines by strong complexes. The region within 8 Mpc from the quasar has been excluded 
because of the proximity effect. No velocity splittings Av < 25 km s -1 were included because 
of the intrinsic line blending due to the typical widths of the Ly a lines. The la deviation from 
a random distribution has been estimated as N~JJ 2 , a good approximation in the case of weak 
clustering £ < 1. The resulting correlation function of the Ly a lines shows a weak but significant 
signal for all three QSOs on scales Av ~ 100 km/s. Exploring the variations of the correlation 
as a function of the column density, an increase appears as the column density threshold is 
raised. The maximum significant signal is obtained for lines with log Njji ^ 13.8, for which the 
measured correlation functions are £2126 = 0.83 ± 0.34 (log Njji > 13.85) £0055 = 0.79 ± 0.24 
(log JVff/ > 13.85), £0014 = 0.81 ± 0.41 (log N HI > 14.1), for At; = 100 km s" 1 (Fig. 1). When 
computing £ for lines in the range log Njji < 13.6 (including in the simulated spectra gaps due 
to lines with log Njji > 13.8), no significant clustering is observed. 

The correlation found for the Ly a cloud positions is less pronounced than for metal-line 
absorption systems or galaxies, but consistent with a scenario of gravitationally induced cor- 
relations, as expected in models where gravitation is an important confining agent ([8]). The 
correlation found would imply, in the standard CDM framework, that the Ly a clouds respon- 
sible for the weakest absorption are abundant in underdense regions. 

Significant anticlustering is detected at scales Av ~ 600 km s -1 (Fig.l). Possible origins 
for this phenomenon are inhomogeneities in the UV ionizing background and variations in the 
baryon to dark matter ratio due to reionization processes (see Bouchet, this meeting). 



2 The dimensions of the Ly a clouds 

At small velocity differences the three-dimensional spatial autocorrelation function £ r and the 
redshift autocorrelation function are related by the convolution ([7]): 

r°° r 00 Hdr f r \ ~ 7 r (Hr-v) 2 (Hr+v) 2 ~ 

( v = Hdr ((r) P(v\r)cx — (-) e"^^ + e" 1 "^ (2) 
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where a Gaussian distribution of the peculiar motions with respect to the Hubble flow is assumed 
together with a power-law spatial correlation function of the galaxy-type. 

At small velocity splittings, the redshift correlation scale r Q depends mainly on the cloud 
sizes r c i and on the velocity dispersion. Although these quantities are poorly known, it is 
instructive to derive constraints on the cloud sizes and velocity dispersions from the observed 
correlation. Assuming an index 7 = 1.77 for the power-law spatial correlation function, the fit 
shown in Fig.l provides upper limits on the cloud sizes as a function of the velocity dispersion 
and of the appropriate correlation scale. We obtain r c i ~ 100 — 190 kpc for a v = 50 — 160 km 
s -1 , and r Q ~ 245 — 420 kpc, respectively. In the former case a good fit to the overall shape of 
the observed correlation function requires a very low value 7 = 1.1. For a v > 200 km s -1 the 
redshift correlation function becomes too flat and very steep 7 values are needed (e.g. a v = 300 
km s -1 requires 7 = 4). 



3 Voids in the Ly a forest 



Voids in the Ly a forest also provide a test for models of the large-scale structure and of the 
homogeneity of the UV ionizing flux. Searches for megaparsec-sized voids have produced a 



couple of claims ([3], [4]), but uncertainties in the line statistics strongly influence the probability 
estimate ([9]). High resolution data, less affected by blending effects, are ideal in this respect. 

We have searched for gaps in a sample of Ly a lines of Q0055-269 limited, to ensure statistical 
completeness and accuracy in the deblending, to log Njji > 13.3 and b < 30 km s -1 . Two regions 
devoid of such lines have been found, centered respectively at A ~ 5000,5206 A, with sizes 
Av ~ 2009,2046 km s -1 (i.e. ~ 20 Mpc). To establish the random probability of observing 
such gaps, a set of 50000 spectra have been simulated with the observed number (N = 178) of 
redshifts randomly generated in the same interval of the data (excluding the region within 8 
Mpc from the quasar) according to the cosmological distribution. The probability of getting a 
gap larger than the maximum observed one is then ~ 0.02. The 2046 km s -1 gap is significant 
only to 2.2cr level, but the joint probability for the presence of the two observed gaps in the 
same spectrum is ~ 2 X 10~ 4 . Even if underdense regions are statistically significant in our 
spectrum, the filling factor is less than 10%. 
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